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We report an analysis of the impressive new lattice simulation results for octet baryon masses 
in 2+1-flavor QCD. The analysis is based on a low order expansion about the chiral SU(3) limit 
in which the symmetry breaking arises from terms linear in the quark masses plus the variation 
of the Goldstone boson masses in the leading chiral loops. The baryon masses evaluated at the 
physical light quark masses are in remarkable agreement with the experimental values, with a 
model dependence considerably smaller than the rather small statistical uncertainty. From the mass 
formulae one can evaluate the sigma commutators for all octet baryons. This yields an accurate value 
for the pion-nucleon sigma commutator. It also yields the first determination of the strangeness 
sigma term based on 2+1-flavor lattice QCD and, in general, the sigma commutators provide a 
resolution to the difficult issue of fine-tuning the strange quark mass. 

PACS numbers: 



In recent years lattice QCD has matured to a level 
where it can be used as a precision tool to confront ex- 
perimental aspects of nonperturbative QCD. This ad- 
vance has occured primarily in the heavy-meson sec- 
tor pp, while for light-quarks, and especially for light- 
quark baryons, progress has been more steady. How- 
ever, new high-precision studies of the baryon spectrum 
in 2+1-fiavour dynamical simulations have recently been 
reported by LHPC |2j, PACS-CS |3], HSC 0] and Durr et 
al. [5]. In this Letter, we demonstrate that these results 
are sufficiently close to the chiral regime that they per- 
mit a reliable determination of the octet baryon masses 
at the physical point, as well the associated sigma terms. 

Our analysis focusses on the (chronologically) first two 
of the above cited studies, those of LHPC and PACS- 
CS. The smallest pion masses are of order 0.29 GeV and 
0.30 GeV, respectively. The lattice spacing for the MILC 
configurations, used by LHPC, has been carefully deter- 
mined in heavy-quark systems [BJ. That same analysis 
also produced a physical determination of the Sommer 
scale, tq = 0.465 ± 0.012 fm, which we use to set the scale 
in the PACS-CS simulations. 

Our working hypothesis is that the (different) im- 
proved actions employed by the two groups yield a very 
good approximation to the continuum theory, so that dis- 
cretization artefacts are small. A unified treatment of the 
continuum extrapolation is only possible once there are 
results available at multiple lattice spacings, such as the 
study of Ref. |S]. On the other hand, a comparison of 
the absolute values of the baryon masses extracted from 
the analysis of the separate LHPC and PACS-CS data 
sets provides a test, a posteriori, of the validity of this 
hypothesis. We shall see that it appears to be a very 
good approximation. 

The first step in our analysis is to calculate the finite- 
volume corrections within the low-energy effective field 
theory (EFT), at each combination of quark masses. 
We can then focus our attention on the extrapolation 



within an infinite- volume, continuum effective field the- 
ory framework. A detailed volume dependence analysis 
of the nucleon mass in 2-flavour simulations has observed 
that the finite-volume corrections are well-described by 
the leading one-loop results of chiral effective field the- 
ory [7]. We extend this approach to the SU(3) case in 
order to estimate the infinite- volume limit of the present 
lattice results. While the leading infrared effects of the 
loop integral are independent of the ultraviolet regular- 
isation, to be conservative we include an uncertainty in 
this correction which amounts to the difference between 
the finite-volume correction determined without a reg- 
ulator [8] and that evaluated with a dipole of mass of 
0.8 GeV [5]. The largest shift we use in our analysis is for 
the nucleon at the lightest pion mass (where m^L = 3.7), 
giving a correction of —0.019 ± 0.005 GeV. 

The chiral expansion of the octet baryons has been 
presented on numerous occasions in the literature — eg. , 
Refs. [HBinilllEGI]- It may be expressed as 

M B = M <°> + 5M B 1] + SM^ /2) + ■■■, (1) 

where the superscript denotes the order of the expansion 
in powers of the quark mass — the explicit chiral sym- 
metry breaking parameter of QCD. Here, denotes 
the baryon mass in the SU(3) chiral limit. 

The SMg^ term is linear in the scale-dependent, run- 
ning quark masses, m; and m s , we refer to Ref. [13] for ex- 
plicit expressions. In this analysis, we replace the quark 
masses by the corresponding meson masses squared, as 
also done in Ref. [2]. To the order that we work in 
this manuscript, the use of either quark masses or me- 
son masses squared is equivalent. 

The next higher order term, 6M B , involves the lead- 
ing one-loop corrections arising from the baryons cou- 
pling to the pseudo-Goldstone bosons, <j> = 7r, K, rj. For 
the baryons, we include both the octet and decuplct. 
Certainly, in the limit -C 6 (the octet-decuplet mass- 
splitting in the chiral limit), the decuplet contributions 
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are formally of higher order. However, with a physical 
kaon mass ~ 0.7 GeV, it is clear that S ~ 0.3 GeV cannot 
be treated as a heavy energy scale. For this reason our 
formal assessment of the order of a given diagram treats 
the octet and decuplet baryons as degenerate. However, 
in order to more accurately represent the branch struc- 
ture in the transition region ~ 5, in the numerical 
evaluation of the loop integrals we maintain the octet- 
decuplet mass splitting. Further, the renormalisation is 
performed such that the decuplet is formally integrated 
out in the limit <C 6. 

For the explicit forms of the loop integrals, we refer the 
reader to Ref . [14] . The coefficients of these loop contri- 
butions [T3] are expressed in terms of the pseudoscalar 
decay constant and the relevant baryon axial charges. 
The nucleoli axial charge, g A = D + F = 1.27 is fixed by 
experiment, while all other couplings are determined by 
SU(6) relations (F = §D and C = -2D). We note that 
C is also related to the decay width of the A, from which 
a similar value can be inferred. We also adopt the chiral 
perturbation theory estimate for the meson decay con- 
stant in the SU(3) chiral limit, / = 0.0871 GeV [15]. The 
octet-decuplet splitting is chosen phenomenologically to 
be the physical N-A splitting, 5 = 0.292 GeV. In prin- 
ciple, all of these input parameters could be constrained 
by actual lattice simulation results, at least in the near 
future — we refer to recent progress in the computation 
of the axial charges [16]. For this study we take these 
parameters from phenomenology, we incorporate gener- 
ous uncertainties in these inputs into our systematic error 
analysis, allowing / to each vary by ±5%, and D, F, C 
and S to vary by ±15%. We leave a global analysis, that 
determines all these inputs simultaneously from the same 
lattice calculation to future work. 

In fitting the lattice results, we wish to minimize the 
uncertainty from higher-order terms in the chiral ex- 
pansion. We therefore limit ourselves to the smallest 
domain of the light-quarks possible with these latest 
lattice results. For both the LHPC and PACS-CS re- 
sults, we include only the results of the simulations for 
ml < 0.2 GeV 2 . 

We first fit the baryon masses to leading- order, keep- 
ing only SM^ in the expansion (llj). This gives a reason- 
able description of the (finite-volume corrected) lattice 
results, with a reduced x 2 (Xdof) °f 0.9 and 0.3 for the 
LHPC and PACS-CS results, respectively. We caution 
reading too much into the absolute value of this naive 
X 2 , as there are certainly strong correlations among the 
points on the same gauge configurations, as mentioned 
in Ref. [2\. In particular, the mass-splittings among the 
baryons will typically be known better than one would es- 
timate from naively adding the uncertainties in the abso- 
lute masses. Without further information on these corre- 
lations, our x 2 wil typically underestimate of the results 
constructed from the full information contained in the 
lattice simulation. 

We now investigate the inclusion of the loop correc- 
tions. We perform fits where the regularisation scale- 



dependence has been removed to all orders and just the 
leading term of the loop integration is retained. Using 
the phenomenological coefficients, the best fit produces 
a (naive) xlof °f order the order 40 (36) for the LHPC 
(PACS-CS) results. Alternatively, a similar fit was also 
done by the LHPC [2], where a suitable Xdof was achcived 
by allowing the chiral coefficients to be fit. While this 
produced a reasonable description of the data, the de- 
termined coefficients (such as g\ or C 2 ) were found to 
be an order of magnitude smaller than phenomenological 
estimates. The most natural conclusion of these obser- 
vations (as also noted in [5]) is that the meson masses 
lie beyond the model-independent or "power counting" 
region (PCR) of the expansion at this order. It is there- 
fore a challenge to maintain the constraints of both the 
EFT and the lattice results (without abandoning one of 
them). 

The breakdown of the expansion at this order should 
come as no surprise, as it has long been known that the 
SU(3) chiral expansion in the baryon sector is quite poor 
- even at the physical quark masses (see eg. [TO] ITT]L 
Donoghuc, Holstcin and Borasoy formulated a solu- 
tion to this problem through long-distance regularisation 
(LDR) [IT]. Concurrently, the same techniques were be- 
ing developed to alleviate the problem of chiral extrapo- 
lation at moderate quark masses in lattice QCD with the 
development of finite-range reg ularisation (FRR) (17) . 
While this has largely been developed for applications 
in SU(2) chiral extrapolations [TB], here we utilise the 
features of the LDR/FRR formalism to perform extrap- 
olations in the framework of chiral SU(3). 

For the relevant formulae describing the renormalisa- 
tion of the loop integrals, we direct the interested reader 
to Ref. JT3]. Upon renormalisation, and to the order we 
are working, the FRR forms produce exactly the same ex- 
pansion as we have described above. The difference lies 
in a resummation of higher-order terms, which are sup- 
pressed by inverse powers of the regularisation scale. At 
this order, the lowest-order addition to the renormalised 
expansion will explicitly appear in the form mi/A, with 
A the regulator mass. 

In working towards ab initio studies of QCD, where no 
external, phenomenological input is used, we do not wish 
to impose any external constraints on the regularisation 
scale. Based on the success of FRR, we should like to 
utilise the observation that the induced resummation of 
higher-order terms provides an improved description of 
a range of lattice obesrvables. Here we use the lattice 
results themselves to determine the regularisation scale 
by minimizing the \ 2 ■ 

Some significant features of using the x 2 measure, with 
the regularisation scale treated as a free parameter, are: 

• no bias of a preferred scale is dictated by phe- 
nomenology; 

• if one is working with results that genuinely lie 
within the PCR, then the x 2 function will be es- 
sentially independent of A; 



3 



• this method provides a quantitative assessment of 
the potential size of the higher-order terms. Fur- 
thermore, an uncertainty arising from the trunca- 
tion of the expansion is automatically incorporated 
into the uncertainties of the fit; 

• if lattice results are included from outside the PCR, 
then through the x 2 1 A is optimised so as to give 
a best-estimate of a resummation of a subset of 
higher-order terms from beyond the working order 
of the chiral expansion. 

To begin, we introduce a single new parameter through 
the regularisation scale. The best-fit to the LHPC 
(PACS-CS) results, with 5 fit parameters, is shown in 
the upper (lower) panel of Fig. [I] where the dipole form 
is shown as an example. As explained above, the fits 
include only those simulation points for mt < 0.2 GeV 2 
(the largest kaon mass is m 2 K ~ 0.40 GeV 2 ). Neverthe- 
less, we see that the level of agreement with the lattice 
simulations at higher m 2 is remarkably good. Further, 
the low-energy constants obtained in the two data sets 
are in agreement to better than half the statistical preci- 
sion of their determination, where for example, the renor- 
malised value, in the SU(3) chiral limit is found to 

be 0.82 ± 0.06 and 0.83 ± 0.08 GeV for the LHPC and 
PACS-CS results respectively. 

The fits are determined by evaluating the baryon mass 
function at each of the lattice kaon and pion masses. The 
curves in Figure [l] are shown for illustrative purposes, 
where the kaon mass at any point on the curve is deter- 
mined by fitting m 2 K as a linear function in m 2 for the 
corresponding lattice ensemble. The physical masses are 
determined by evaluating the fit function at the physical 
pion and kaon masses (which has no bearing on the linear 
form used for the figure). To illustrate how this extrap- 
olation in the strange-quark mass works, the lower panel 
of Fig. [I] shows a fit to just the two PACS-CS results at 
fixed k s . Using that fit, the results of the simulation at 
the different k s are shown as a prediction — by evaluating 
the fit function at the lattice kaon mass of this ensemble. 
The agreement between the results of the simulation and 
the predicted values is illustrative of the reliability of the 
fit in estimating the dependence of the octet masses on 
the strange-quark mass. In our final results, the lattice 
points at this extra strange-quark mass are also included. 

For the LHPC (PACS-CS) results the optimal dipole 
regularisation scale is found to be A = 1.1 ± 0.4 GeV 
(0.91 ± 0.34 GeV). The minimum xLt is °- 25 (0.05), 
where the improvement over the above linear forms is 
evident. As discussed, the existence of a preferred regu- 
larisation scheme is a direct signature that the results lie 
outside the PCR. Nevertheless, the extrapolated preci- 
sion obtained is rather encouraging, considering that the 
uncertainties incorporate the effect of the relatively large 
range permitted for the regularisation scale, as much as 
~ 0.7-1.5 GeV. 

Given that our results do demonstrate that we are out- 
side the PCR, we also investigate the model-dependence 
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FIG. 1: The upper and lower panels show, respectively, the 
dipole fits (lines) to the LHPC and PACS-CS lattice simula- 
tion results (circles) of the octet baryon masses (curves top 
to bottom show S, S, A and N). The squares at the physi- 
cal pion mass display the extrapolation to the physical quark 
masses. The stars denote the physical baryon masses. The 
errors indicated by the bands represent the total statistical 
errors, including the variation of A. The third lightest pion 
mass of the PACS-CS results, which involved a lower strange 
quark mass, is not included in the displayed fit. This allows 
a prediction (squares) based on the other points (see text). 



in the choice of regulator. Smooth forms, such as a 
dipole, monopole and a Gaussian give indistinguishable 
results, with the most significant difference to these seen 
in a sharp cutoff, as was seen in Ref. [13]. Our esti- 
mate of the model-dependent uncertainties incorporates 
the full variation over these different functional forms, in 
addition to all the uncertainties in the phenomenological 
input parameters, as described above. 

In Table |T] we report best-estimates of the physical 
masses by combining the different lattice simulations as- 
suming they are each good approximations of the con- 
tinuum limit. We treat the difference between the sim- 
ulations as an estimate of the discretisation uncertainty. 
The size of these estimates are consistent with the find- 
ings of Ref. [5] . We also report the dimensionless baryon 
sigma terms, a Bq = (m q /M B )dM B /dm q . We note that 
the model-dependence, is relatively small compared with 
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TABLE I: Extracted masses and sigma terms for the physical 
baryons. The first error is statistical; the second estimates the 
discretisation artifacts by the difference between the results 
for the LHPC and PACS-CS results; the third error represents 
a model-dependence estimate as described in the text. A fur- 
ther ~ 2% uncertainty in the absolute scale [6] is understood. 
The experimental masses are shown for comparison. 



Mass (GeV) Expt. 



0.945(24) (4) (3) 0.939 
1.103(13) (9) (3) 1.116 
1.182(11) (2) (6) 1.193 
1.301(12)(9)(1) 1.318 



0.050(9) (1) (3) 0.033(16) (4) (2) 

0.028(4) (1) (2) 0.144(15) (10) (2) 

0.0212(27)(1)(17) 0.187(15)(3)(4) 

0.0100(10)(0)(4) 0.244(15)(12)(2) 



the present statistical precision. This model-dependence 
will not decrease with increased statistics and therefore 
the present analysis, at this order and at these quark 
masses, is precision limited. Improvements beyond this 
will certainly require increased efforts in the EFT and 
numerical computations deeper in the chiral regime. 

Our results show agreement between the absolute val- 
ues of the baryon masses and the corresponding experi- 
mental values. This is an important independent confir- 
mation of the absolute scale determination reported by 
Aubin et al. [5]. Further, the sigma terms extracted by 
differentiating the fitting formulae, have some quite in- 
teresting features. The pion-nucleon sigma commutator 
is, within uncertainties, compatible with phcnomcnologi- 
cal estimates [H]. The strangeness sigma commutator is 
consistent with best EFT estimates [10], yet an order of 
magnitude more precise. This small value is observed to 
be consistent with "unquenching" estimates , as well 
as a recent 2-flavour lattice QCD estimate [2"T] . 

In line with quark-model expectations, we note that 
the 2 is most sensitive to the strange quark mass. The 
linear projections of kaon masses of the LHPC and PACS- 



CS results, used in Fig. [TJ give strange quark masses (or 
rn %~ m ir /2) that are ~ 30% too high. With the definition 
of the sigma term giving = (JBq-^, and as s ~ 0.24 
from Table [TJ this indicates that the 3 mass on each curve 
should lie high by roughly 0.24 x 0.30 ~ 7%. 

In summary, we have demonstrated an excellent de- 
scription of current 2+1-flavour lattice QCD results 
based on a low-order SU(3) chiral expansion. While the 
expansion to this order is not sufficiently convergent to 
ensure that the fits are model-independent, we find that 
the model-dependence is actually small compared with 
the current statistical precision. In the future, as more 
simulations and statistics accumulate, the significance of 
the model-dependent component of the error will increase 
and will necessitate an increased effort in the EFT. We 
also anticipate that future studies will strengthen the nu- 
merical constraints on both the continuum and infinite- 
volume extrapolations, such as the work of Rcf . |5j. 

Within the caveats of the present study, we have 
demonstrated a robust and precise determination of the 
absolute values of the octet bayon masses. The controlled 
extrapolations have also permitted a reliable determina- 
tion of the baryon sigma terms, where we have been able 
to extract an accurate determination of the a commuta- 
tor as well as showing that the strange-quark sigma term 
of the nuclcon is considerably smaller than phcnomcno- 
logical estimates. As just one example of the importance 
of these results, we note the significance for dark matter 
searches explained in Ref. [22 . 
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